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ABSTRACT 

We report the detection of a 0.6 Mj extrasolar planet by WASP-South, WASP-25b, 
transiting its solar-type host star every 3.76 days. A simultaneous analysis of the 
WASP, FTS and Euler photometry and CORALIE spectroscopy yields a planet of 
Rp = 1.22 Rj and M p = 0.58 Mj around a slightly metal-poor solar-type host star, 
[Fe/H] = -0.05±0.10, of R* = 0.92 R Q and M* = 1.00 M . WASP-25b is found to have 
a density of p p = 0.32 pj, a low value for a sub- Jupiter mass planet. We investigate 
the relationship of planetary radius to planetary equilibrium temperature and host 
star metallicity for transiting exoplanets with a similar mass to WASP-25b, finding 
that these two parameters explain the radii of most low-mass planets well. 
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1 INTRODUCTION 

To date, over 440 exoplanets have been discovered, including 
more than 70 detected by the transit methoo3. The tran- 
sit method together with follow-up radial velocity observa- 
tions allow measurement of both the mass and radius of 
the planet, leading to a v alue for the planet's bulk density 
jCharbonneau et aljlioool ). The atmospheric composition of 
transiting exoplanets can also be investigated through high- 
pre cision photometric and spe ctroscopic measurements, see 
e.g. ICharbonneau et all (|2002t ). 

A wide range of transiting exoplanets radii has been 
found and there has been much investigation into the fac- 
tors that may infl uence a planet's radius. For example, 
iGuillot et "all l|2006l ) propose a negative relationship between 
the metallicity of a host star and the radius of an orbit- 
ing planet, caused by an increase in the amount of heavy 
elements in the planet, leading to a more massive core 
and hence smaller radius for a given mass. Alternatively, 
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iBurrows et all (|2007l ) consider that increasing the metallic- 
ity may increase the opacity of an exoplanet's atmosphere, 
retarding cooling and leading to a larger radius for a given 
mass. Another influence on a planet' s radius may be the 
equil ibrium temperature of the planet (|Guillot fc Showman! 
I2002T ). determined by the stellar irradiation and the planet's 
distance from its host star. Tidal heating due to the circu- 
larisation of the orbits of close-in exopl anets may also play 
a rol e in inflating the plan etary radius (|Bodenheimer et al.l 
120031 ; IJackson et alj I2OO8I ). One motivation of the Super- 
WASP project is to detect enough transiting exoplanets, 
with a wide range of orbital and compositional parameters, 
to allow analyses that may distinguish between such differ- 
ing models. 

In this paper, we report the discovery of a 0.6 Mj 
planet orbiting a solar-mass star, WASP-25 (=TYC6706- 
861-1, =1SWASP J130126.36-273120.0), in the southern 
hemisphere. Analysis of photometric and spectroscopic data 
reveals WASP-25b to b e another low-density pla net, com- 
parable to HD 209458b (|Charbonneau et al.lf2000h . We also 
analyse the dependence of the radii of low-mass planets on 
host star metallicity and planetary equilibrium temperature, 
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Figure 1. WASP discovery lightcurvc folded on the orbital period 
of P = 3.765 d. Points with error above 3 X median were clipped, 
where median = 0.012 mag 

including WASP-25b and 18 other transiting planets, finding 
a relationship using Singular Value Decomposition analysis 
that gives an excellent agreement between observed and cal- 
ibrated radii. 

In Section [2] we describe the photometric and spec- 
trscopic observations and data reduction procedures. In Sec- 
tion [3] we present the stellar and planetary parameters ex- 
tracted from these data. Finally, in Section [4] we compare 
WASP-25b with the ensemble of known planets of similar 
mass, and examine the relationship between stellar metal- 
licity, irradiating flux (via the planet's equilibrium temper- 
ature) and planet radius. 
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Figure 2. Top lightcurve: FTS follow-up photometry of WASP- 
25 on 3 April 2010. The central transit time is at HJD = 
5290.05617. Bottom lightcurve: Euler photometry on 10 April 
2010 with central transit time at HJD = 5297.58552. 



2 OBSERVATIONS 

2.1 Photometric Obervations 

The WASP-South observatory is located at SAAO in South 
Africa, and consists of eight 11cm telescopes of 7.8° x 7.8° 
field of view each, on a single fork mount. The cameras scan 
repeatedly throug h eight to ten sets of fi elds, taking 30 sec- 
ond exposures. See lPollacco et al.l (|200rJ ) for more details on 
the WASP project and the d ata reduction procedure, and 
ICollier Cameron et al.l (|2007h and lPollacco et al.l (|2008l ) for 
an explanation of the candidate selection process. 

WASP-25 was observed by WASP-South in 2006, 2007 
and 2008, producing a total of 14,186 photometric data- 
points. The 2007 dataset showed the transit event most 
clearly, detected at a period of 3.76 days. Figure [T] shows 
the WASP-South discovery lightcurve, using data from all 
seasons. 

Further photometric observations were subsequently 
obtained on 3 April 2010 using the LCOGT 2m Faulkes Tele- 
scope South (FTS) at Siding Spring, Australia. Observations 
were obtained using the fsOl Spectral camera containing a 
4096 x 4096 pixel Fairchild CCD which was binned 2x2 giv- 
ing 0.303 arcsec/pixel and a field of view of ~ 10' x 10'. 285 
datapoints were obtained through a Pan-STARRS z filter, 
capturing an entire transit. 

The data were pre-pro cessed through the WASP 
Pipeline (|Pollacco et al.l [20061 ') to perform masterbias/flat 
creation, debiassing and flatfield correction in the standard 
manner. Object detection and ape rture photom etry were 
performed using the DAOPHOT (|Stetsonl 1 19871 ) package 



within the IRAF environment with an aperture size of 
18 binned pixels. Differential magnitudes of WASP-25 were 
formed by a weighted combination of the flux relative to 26 
comparison stars within the field of view 

A second full transit was observed a week later on 10 
April 2010 in the R-band with EulerCam on the 1.2m tele- 
scope at La Silla, Chile, comprising 204 datapoints. The 
resulting lightcurves are shown in Figure [2] 



2.2 Spectroscopic Observations 



WASP-25, a V„ 



11.9 star, was observed 29 times with 



the CORALIE spectrograph on the 1.2m Euler telescope, 
between 29 December 2008 and 28 June 2009. The spectra 
were processed using the standard data reduction pi peline 
for CORALIE (|Baranne et al.lll996l ; lMavor et al.ll2009h. plu s 
a correction for the blaze function (jTriaud fc et a] H2010T l. 
These data are given in Table [T] and shown phase- folded in 
Figure [3] the low-amplitude radial velocity variations and 
the lack of correlation betwe en the bisector spans and radial 
velocity jQueloz et alj|200lh . shown in Figure [4] are consis- 
tent with a planet-mass object orbiting the host star. 



2 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Table 1. CORALIE radial velocity measurements of WASP-25 



BJD-2 400 000 


RV 


°RV 


BS 




(km s- 1 ) 


(km s- 1 ) 


(km s _1 


54829.8227 


-2.577 


0.013 


-0.019 


54896.7698 


-2.651 


0.011 


-0.022 


54940.7092 


-2.716 


0.012 


0.010 


54941.7043 


-2.619 


0.012 


-0.012 


54942.7257 


-2.576 


0.012 


-0.029 


54943.6374 


-2.618 


0.012 


-0.055 


54944.7155 


-2.680 


0.012 


-0.032 


54945.7265 


-2.615 


0.013 


-0.011 


54946.6166 


-2.582 


0.013 


-0.019 


54947.6016 


-2.641 


0.011 


-0.000 


54947.7912 


-2.689 


0.013 


0.001 


54948.6130 


-2.704 


0.011 


-0.011 


54949.8031 


-2.551 


0.018 


-0.030 


54950.6221 


-2.591 


0.013 


-0.044 


54951.6953 


-2.701 


0.012 


-0.011 


54971.6453 


-2.678 


0.021 


0.026 


54972.6724 


-2.561 


0.013 


0.027 


54973.5157 


-2.587 


0.013 


-0.030 


54974.6787 


-2.714 


0.014 


-0.064 


54975.5379 


-2.667 


0.014 


-0.023 


54976.6837 


-2.556 


0.013 


-0.011 


54982.6194 


-2.664 


0.021 


-0.029 


54983.6213 


-2.568 


0.015 


0.043 


54983.6446 


-2.597 


0.015 


-0.029 


54984.5785 


-2.558 


0.015 


0.023 


54985.6100 


-2.699 


0.012 


0.004 


54995.5555 


-2.509 


0.014 


-0.061 


55009.6287 


-2.606 


0.018 


0.018 


55010.5967 


-2.539 


0.023 


-0.036 
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Figure 3. Top plot: Radial velocity measurements. The solid line 
is the best-fitting MCMC solution. The centre-of-mass velocity, 7 
= -2.6323 km s — 1 , was subtracted. 



3 SYSTEM PARAMETERS 
3.1 Stellar Parameters 

The CORALIE spectra were placed on a common wave- 
length scale and co-added to produce a higher signal-to-noise 
spectrum, allowing an analysis of the host star and thus mea- 
surement of the stellar temperature, gravity, metallicity, v 
sin i and elemental abundances, given in Table [2] where 
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Figure 4. Bisector spans versus radial velocity, where bisector 
uncertainties are taken to be equal to twice the radial velocity 
uncertainties. 
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Table 2. Stellar Parameters of WASP-25 



Parameter 


Value 


T cff 


5750 ± 100 K 


log g 


4.5 ± 0.15 


V 


1.1 ± 0.1 km s" 1 


v sin i 


3.0 ± 1.0 km s" 1 


Spectral Type 


G4 


[Fe/H] 


-0.05 ± 0.10 


[Si/H] 


0.00 ± 0.06 


[Ca/H] 


0.08 ± 0.14 


[Ti/H] 


0.04 ± 0.07 


[Ni/H] 


-0.08 ± 0.10 


log[Li/H] 


1.63 ± 0.09 


V (mag) 


11.88 



SWASP J130126.36-273120.0 



77 is microturbulence. The spectral data wer e analysed us- 
ing t he UCLSYN spectral synthesis package |Smallev et al.l 
1200 ll) and ATLAS9 models without convective overshooting 
( Castelli et ai1ll997f ). The effective temperature and log g 
were determined using the Ha and H/3 lines, and the Na ID 
and Mg lb lines respectively, while the Ca H and K lines pro- 
vided a check on th ose values. Further details of the spectral 
analysis are given in I West et al.l (|2009h - The analysis yielded 
T eff = 5750 ± 100K and [ Fe / H l = -0.05±0.10. 



3.2 System Parameters 

The WASP-South, FTS and Euler photometry were simulta- 
neously analysed with the CORALIE radial velocity data in 
a Markov-Chain Monte-Carlo (MCMC) analysi s. This anal- 
ysis is described in lCollier Cameron" et all (120071 ) . but is here 
modified as described in lEnoch et all i|2010t ) to determine 
stellar mass using a calibration on T g g-, log p and [Fe/H] 
(si milar to the T„ g , log g and [Fe/H] calibration described 
in iTorres et all (12111)1)1)). The temperature and metallicity 
were obtained through spectral analysis (given in Table [2}, 
and the density is determined directly from the photometry. 
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Figure 5. Isochronc tracks from Girardi et al. (2000) for WASP- 
25 using the best-fit mctallicity of -0.07 and stellar density 1.29/9© . 
Isochrones are: 0.01995 (solid), 0.7943 (dashed), 1.99 (dot-dash), 
3.16 (dotted), 5.01 (triple dot-dash) Gyr. 



Mass and radius values obtained for other WASP host stars 
via this method agree closely with values obtained through 
isochrone analysis. 

An initial analysis was performed allowing the eccentric- 
ity value to float, producing a value of e = 0.123lQ Q4g. This 
eccentricity value is consistent with at the 3<r level, and was 
suspected to not be significant, as can occu r from spurious 
asym metries in quadrature fits due to noise (|Laughlin et al.l 
l2005l ), so a second analysis was performed with eccentricity 
fixed to 0. We performed an F-test on the radial velocity 
residuals from the circular and floating eccentricity fits, re- 
sulting in a value of 0.897 which shows that the eccentric fit 
is not significant. 

The resulting circular model best-fit parameters for the 
star-planet system are listed in Table O using the best-fit 
parameters from the e = fit and the uncertainties from 
the fit allowing eccentricity to float, to sufficiently account 
for uncertainty due to unknown eccentricity. The results 
show WASP-25 to be a solar analogue of one solar mass and 
0.92 ± 0.04 solar radius, and WASP-25b to be a bloated 
hot Jupiter of 0.58 ± 0.04 Mj and 1.22 ± 0.06 Rj, giving 
a planet density of p = 0.32^o'q3 pj. 

The stellar density of 1.29±— 0.11 p® obtained from the 
MCMC analysis was used along with the determined stel- 
lar te mperature and meta llicity values in an interpolation 
of the iGirardi et all (|2000T ) stellar evolution tracks, see Fig- 
ure [5] Using the best-fit metallicity of —0.07 indicates that 
WASP-25 has a mass of 1.02±q ™M q , agreeing well with 
the calibrated MCMC result of 1.00 ± O.O3M , and an age 
of 0.02±g;g? Gyr. 
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Figure 6. Radius versus metallicity of 18 planets of 0.1-0.6 Mj. 
4 DISCUSSION 

A density of 0.32l[5:o3 PJ places WASP-25b amongst the 
bloated hot Jupiters, with ov er 80% of known tra nsiting ex- 
oplanets being more dens<B iGuillot et all (|2006l ) proposed 
a correlation between the metallicity of a host star and the 
amount of heavy elements present in the planet. A larger 
amount of heavy elements is likely to produce a more mas- 
sive rocky core and hence would lead to larger planetary 
radii for lower met allicity stars, for a given planetary mass 
(|Fressin et al.|[2007h . Alternatively, the heavy elements could 
increase the opacity of the planetary inte rior, potentially in- 
creasing the radius l|Burrows et al.ll2007l ). 

WASP-25b has a mass close to the value above which 
the gravitational, potential energy begins to exceed the elec- 
trical binding energy of electrons in atoms as the main force 
opposing electron degeneracy pressure i n the deep planetary 
interior (|Lvnden-Bell fc O'Dwyeil 120011 ). It is in this lower 
mass range that we see the greatest range of radii at a given 
mass. This suggests that environmental and compositional 
differences between individual planets have the greatest in- 
fluence on their radii in this mass range. 

To investigate this relationship, we calculated the corre- 
lation between radius and metallicity for all known transit- 
ing exoplanets for which metallicity data was available (67 
planets) , finding only a very weak relationship with a corre- 
lation coefficient of -0.17. However, splitting the planets into 
low-mass (0.1 — 0.6 Mj) and higher-mass (> 0.6 Mj) sets re- 
sulted in correlation coefficients of -0.53 (with a probability 
of no correlation of 0.022) and 0.07 respectively, indicating 
that metallicity plays a role in determining the radii of low- 
mass planets but not those of higher mass. We plot the radii 
against host star metallicity of low-mass planets in Figure |S1 
We concentrate on the set of 18 low-mass planets, of which 
9 are WASP planets, including WASP-25b, see Table gj), in 
further investigating planetary radii. 

A major influence on the radius of an exo planet is likely 
to be the planet's equilibrium temperature l|Fressin et al.l 
I2007D , T eq , defined as 

3 www.exoplanet.eu 
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Table 3. System Parameters of WASP-25 



Parameter 



Symbol 



Value 



Period (days) 
Transit Epoch (HJD) 
Transit duration (days) 
Planet/Star area ratio 
Impact Parameter 



P 3.764825 ± 0.000005 

T 5274.99649 ± 0.00017 

D 0.116 ± 0.001 

Rl/Rl 0.0187 ± 0.0002 

b 0.38+»'°« 



Stellar Reflex Velocity (ms -1 ) 
Centre-of-mass Velocity (ms -1 ) 
Orbital separation (AU) 
Orbital inclination (deg) 
Orbital eccentricity 



ifi 75.5 ± 5.3 

7 -2632.3 ± 0.6 

a 0.0473 ± 0.0004 

i 88.0 ± 0.5 

e (adopted) 



Stellar mass (Mq) 


M» 


1.00 


± 


0.03 


Stellar radius (Rq) 


R* 


0.92 


± 


0.04 


Stellar surface gravity (log go) 


log g* 


4.51 


± 


0.03 


Stellar density (pq) 


P* 


1.29 


± 


0.10 


Stellar metallicity 


[Fe/H] 


-0.07 


± 


0.10 


Stellar effective temperature 


T e// 


5703 


± 


100 



Planet 


mass (Mj) 


M p 


0.58 ± 0.04 


Planet 


radius (Rj) 


Rp 


-, 99 +0.06 
1 " a ' 5 -0.05 


Planet 


surface gravity (log gj) 


log 9p 


2.95 ± 0.04 


Planet 


density (pj) 


Pp 


n oi 7 +0.036 
u.oi / „o.031 


Planet 


temperature (A=0, F=l) (K) 


Teq 


1212 ± 35 



where T e g- is the host star's effective temperature, R* is 
the stellar radius, a is the semi-major axis, A is the planetary 
Bond albedo and F is the fraction of the planet re-radi ating 
flux. Jupiter has a Bond albedo of 0.28 (|Tavlorlll965l ) but 
albedo values have not so far been determ ined for the ma- 
jority of exoplanets, and I Ro we et all (113) found a very low 
albedo for HD 209458b, at just 0.04 ± 0.05. We set A = 
(F = 1) here to calculate the equilibrium temperatures. Not 
knowing the true albedo values for all the planets studied 
here produces an uncertainty in the values calculated for 
equilibrium temperature of (1 — A) 1 ^ 4 . We find the correla- 
tion coefficient between the (approximate) equilibrium tem- 
perature and radius of all 18 planets to be 0.55, a moderate 
relationship, shown in Figure [7] The correlation of radii and 
equilibrium temperature is stronger than this for the major- 
ity of planets - removing the two outliers, WASP-17b and 
HD 149026b, from the correlation produces a coefficient of 
0.82 for the remaining 16 planets, a strong relationship. Ac- 
curate albedo values, and hence equilibrium temperatures, 
may produce an even stronger correlation than found here. 

Since both the host star metallicity and the planet's 
equilibrium temperature seem to explain the radii of low- 
mass planets in part, we performed a weighted Singular 
Value Decomposition (SVD) fit to fully quantify the rela- 
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Figure 7. Radius versus equilibrium temperature of 18 planets 
of 0.1-0.6 Mj. 

tionship. This resulted in the equation 

R p = 0.48 + 0.50(T e q/1000A') - 0.61[Fe/H] (2) 

producing the fit to the planetary radii shown in Figure [8] 
Removing the two outliers mentioned above, WASP-17b and 
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Figure 8. Results of SVD fit on the radii of 18 planets of 0.1- 
0.6 Mj. 
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Figure 9. Results of SVD fit on the radii of planets of 0.1-0.6 Mj, 
excluding WASP-17b and HD 149026b. 

HD 149026b, leads to a much better fit, shown in Figure O 
implying that additional parameters particularly affect these 
two planets. The SVD fit for the 16 planets is 

R p = 0.20 + 0.73(Teq/1000iO - 0.25[Fe/H] (3) 

The negative slope of the metallicity correlation and of 
the SVD coefficient implies that the effect of a massive core 
(|Guillot et al.|[2006l ) outweighs the effect of increased opacity 
in the interior of the planet, i.e. increased metallicity leads 
to a smaller planetary radius. 



5 CONCLUSIONS 

We have reported the detection of a 0.58 Mj planet, WASP- 
25b, transiting a slightly metal-poor solar-mass star in the 
southern hemisphere with an orbital period of 3.76 days. 
WASP-25b has a low density, 0.29 pj, and we investigate 
its bloated radius, R p — 1.26 Rj. We find that the radii 
of most transiting exoplanets of a similar mass to WASP- 
25b can be explained well by a calibration to the host star 
metallicity and planetary equilibrium temperature. 
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